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29 ■ 1 Introduction 

On 

To understand the formation of the visible structure in the Universe out of 
i ■ an initially almost homogeneous matter distribution is one of the most fas- 

cinating quests of modern cosmology. The first step to this understanding is 
| of course an assessment of the matter distribution in the present Universe on 

very large scales extending over several hundred Mpc (for a Hubble constant 
of 50 km s _1 Mpc -1 ). Such large scales are interesting for two major rea- 
sons: the present day structures on this scale are directly comparable to the 
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5_i ■ signature of the primordial structure detected in the microwave background 



by COBE with similar co-moving sizes, and at these scales the observed 
density fluctuations are just linear amplifications of the initial conditions in 
the early Universe. While the study of the galaxy distribution has already 
given us very interesting insights into the structure on scales of a few hun- 
dred Mpc an alternative approach using the next larger building blocks of 
the Universe, galaxy clusters, as probes for the cosmic structure can give us 
a ready access to even larger scales (see also Tadros et al. 1998 for the APM 
cluster survey in the optical). 

X-ray astronomy has offered a unique tool to efficiently detect and char- 
acterize galaxy clusters out to large distances. Originating in the hot intra- 
cluster plasma that fills the gravitational potential well of the clusters, the 
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X-ray emission is an equally robust parameter for a first estimate of the size 
and mass of clusters as the velocity dispersion measured from the galaxy 
redshifts. But, while the X-ray luminosity can be readily detected for many 
clusters in an X-ray sky survey, the collection of velocity dispersions requires 
very large and time consuming redshift surveys. 

In the project described here we have embarked on a redshift survey 
of galaxy clusters detected in the ROSAT All-Sky Survey (Triimper 1993, 
Voges et al. 1996), the first complete sky survey conducted with an imag- 
ing X-ray telescope. In the frame of an ESO key programme (Bohringer 
1994, Guzzo et al. 1995) we are seeking a definite identification of all possi- 
ble cluster candidates found in the ROSAT Survey in the southern celestial 
hemisphere above an X-ray flux limit chosen such to provide a reasonably 
homogeneous sensitivity coverage of the sky area. Within this survey pro- 
gramme, which we call the ROSAT ESO Flux Limited X-ray (REFLEX) 
Cluster Survey, we are investigating about 800 galaxy clusters and have 
obtained more than 300 new cluster redshifts. 

The present article describes methods of the cluster detection and iden- 
tification as well as first preliminary results on the statistics of the cluster 
population and a view on the large-scale distribution of the clusters. In a 
following paper we describe in more detail the observing strategy and the 
spectroscopic results. 

2 Cluster Detection and Follow-Up Observations 

The X-ray sky atlas constructed from the ROSAT All-Sky Survey Mission 
with its more than 100,000 X-ray sources (Voges et al. 1996) contains thou- 
sands of mostly unidentified galaxy clusters. Since the majority of these 
sources are characterized by few photons, only the brightest, well extended 
X-ray cluster sources are readily identified, while the main part of the iden- 
tifications has to be based on further optical information. For the REFLEX 
Survey we made mainly use of the COSMOS optical object catalogue (e.g. 
Heydon-Dumbleton et al. 1989) produced from the digital scans of the pho- 
tographic UK Schmidt Survey ( providing star /galaxy separation for the 
sky objects with high completeness down to bj = 20.5 mag.). The cluster 
candidates are found as overdensities in the galaxy density distribution at 
the position of the X-ray sources (see Bohringer et al. in preparation). Not 
all the sources flagged by this technique are true galaxy clusters, however. 
The price paid for aiming at a high completeness in the final catalogue is 
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Figure 1: One of the galaxy clusters with a redshift of z = 0.188 not listed in 
the catalogue by Abell et al. (1989) found during the identification of X-ray 
cluster sources from the ROSAT All-Sky Survey in the frame of the ESO key 
programme. The figure has been produced from the STScI digital scans of 
the UK-Schmidt IllaJ plates with X-ray cotours from the ROSAT All-Sky 
Survey images superposed. The contours are 1.5,2,3,4,... sigma significance 
contours of the X-ray signal accounting properly for the photon noise and 
backgound. 

a low detection threshold in the galaxy density leading to a contamination 
of the candidate list by more than 30% non-cluster sources. This large 
contamination can be reduced to about one third by a direct inspection of 
the photographic plates, the detailed X-ray properties, and using all the 
available literature information. 

The subsequent work is the observational task of the ESO key pro- 
gramme comprising a total of 90 observing nights distributed evenly among 
the 1.5m, 2.2m, and 3.6m telescopes at La Silla. These follow-up observa- 
tions have a two-fold purpose: we search for a definite identification for the 
unknown objects and collect redshifts for all clusters for which this param- 
eter is not known. To make optimal use of the ESO telescopes we observe 
nearby, poor clusters and groups in single-slit mode with the smaller tele- 
scopes and use the 3.6m telescope with EFOSC in multi-slit operation to get 
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multiple galaxy spectra of dense, rich clusters. Several coincident redshifts 
strongly support the cluster identification and the spectroscopy of central 
dominant early type galaxies at the X-ray maximum plays a particularly 
important role in this study. 

The observing programme will be completed in fall of this year. Having 
observed the clusters with higher X-ray flux with highest priority, we can 
construct a first catalogue of clusters for the brighter part of the sample 
down to an X-ray flux limit of 3 • 10~ 12 erg s _1 cm -2 (in the ROSAT band 
0.1 - 2.4 keV) comprising 475 objects. 53% of these clusters can be found 
in the main catalogues of Abell (1958) and Abell, Corwin, & Olowin (1989) 
and further 10% in the supplementary list. Most of the remaining clusters 
are previously unknown objects. This result highlights the importance of 
the selection process in the construction of the sample: a significant fraction 
of the X-ray bright and massive clusters would have been missed had the 
X-ray clusters only been identified on the basis of existing catalogues. The 
fraction of known Abell clusters decreases further if one goes to lower X- 
ray flux limits e.g. in the extended REFLEX sample. We should also note 
here, however, that one part of the objects missing in Abell's compilation 
are nearby, X-ray bright groups, which are not rich enough in their galaxy 
content to fulfill Abell's criteria. 

The high completeness of this sample is demonstrated by a counter-test, 
in which we searched for X-ray emission at the position of all Abell and 
ACO clusters in the ROSAT Survey independent of a previous detection of 
these sources by the survey source identification algorithm. Only 5 clusters 
(~ 1%) were found to have been missed by the cluster search based on the 
COSMOS data. Fig. 1 gives an example of a non-Abell clusters found at a 
redshift of z = 0.188. A fraction of the X-ray sources in the cluster candidate 
list are found to have non-cluster counterparts in the follow-up observations 
including AGN in clusters where the AGN provide the main source of the 
X-ray emission. 

A few spectacular discoveries were made in the course of this programme 
including the most X-ray luminous cluster discovered so far, RXJ1347-1145 
(Schindler et al. 1995,1996) and clusters featuring bright gravitational arcs, 
e.g. S295 (Edge et al. 1994). 
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Figure 2: Cumulative number counts of X-ray clusters as a function of the 
limiting flux for the REFLEX sample. 

3 Properties of the X-ray Clusters of Galaxies 

In the following we are reporting results based on the X-ray bright sample 
of 475 galaxy clusters for which 413 cluster redshifts have been obtained and 
reduced so far. A plot of the number counts of the cluster population as a 
function of the limiting X-ray flux is shown in Fig. 2. The logarithmic graph 
has a slope that is in perfect agreement with an Euclidian slope of —3/2. 
This slope is easily explained by the fact that the majority of the clusters 
is not very distant - the peak in the redshift histogramme is at z ~ 0.06 - 
and thus band-corrections and evolutionary corrections are not important. 
Fig. 3 shows the distribution in redshift and X-ray luminosity. While most 
of the clusters detected are not very distant (median redshift is z = 0.085), 
a few very luminous clusters are found out to redshifts of z = 0.3 with 
one outstanding object at a redshift of z = 0.45, the most luminous cluster 
mentioned above. One also notes clearly the break at a redshift of z = 0.3, 
which is caused by the limited depth of the optical plate material used for 
the cluster pre-identification. Thus we can expect that there are more, in 
fact very interesting luminous clusters in the X-ray source list for this flux 
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Figure 3: Distribution of the X-ray luminosities as a function of redshift for 
the galaxy clusters of the REFLEX sample. 

limit with redshifts between z = 0.3 to 0.5. But a more expensive imaging 
search programme would be necessary to find these objects. 

Another very interesting part of this X-ray cluster population is found 
among the nearby, low X-ray luminosity objects. These are groups of galax- 
ies dominated by very giant cD galaxies which are often found to be two or 
three magnitudes brighter than the next brightest galaxy in the system. A 
CCD image of such a group is shown in Fig. 4. 

The most straightforwardly obtained and important distribution func- 
tion of the cluster sample is the X-ray luminosity function. This function 
is most closely related to the mass function of the clusters which is used as 
an important calibrator of the amplitude of the density fluctuation power 
spectrum of the Universe (e.g. White et al. 1993). A preliminary version of 
the X-ray luminosity function of the sample is shown in Fig. 5. Note that 
this function was derived for the sample when (~ 80%) of the redshifts had 
been obtained. But in spite of this incompleteness the luminosity function 
already recovers the densities reached in previous surveys (e.g. Ebeling et 
al. 1997, DeGrandi 1996) to which the present result is compared in the 
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Figure 4: Galaxy group with a very dominant central cD galaxy at a redshift 
of z = 0.0313 in the REFLEX sample. The CCD image is a short routine 
exposure (without filter) to define the MOS slit mask taken with EFOSC1 
at the 3.6m telescope. The scale of the image is about 3.5 by 3.5 arcmin. 

figure. The high quality of the present sample is also demonstrated by the 
fact that the slope of the number count function of Fig. 2 is Euclidian while 
the previous surveys showed a significantly lower value for the logarithmic 
slope of ~ —1.35. 

As pointed out in the introduction, the most exciting aim of this pro- 
gramme is to assess the large-scale structure of the galaxy cluster and matter 
distribution in the Universe. The currently most popular and important sta- 
tistical function for the characterization of the large-scale structure is the 
density fluctuation power spectrum. A first preliminary result for this func- 
tion obtained from 188 cluster of our sample in a box of 400 Mpc is shown 
in Fig. 6 (for details see Schuecker et al. in preparation). The function is 
featuring an interesting maximum at a scale of about 100 - 150 /i^qq Mpc. 
The location of this maximum is related to the size of the horizon of the 
Universe at an epoch when the energy density of matter and radiation was 
equal and it is therefore a very important calibration point for the theoret- 
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Figure 5: Differential X-ray luminosity function for the REFLEX sample 
(points with error bars; for details see Bohringer et al. in preparation) 
compared to the previous results by DeGrandi (1996) and Ebeling et al. 
(1997). A value of Hq = 50 km s _1 Mpc -1 is assumed for the scaling. 



ical modeling of the evolution of our Universe. Note that the construction 
of the power spectrum shown in Fig. 6 is based on a luminosity (mass) 
selection of the clusters which varies with redshift. Therefore the quantita- 
tive interpretation of the this result is not straightforward and more work is 
needed to relate this function to the power spectrum of the matter density 
fluctuation in the Universe. 

From the power spectrum shown we can also derive another statistical 
function which is more illustrative to non-experts, the rms-fluctuation level 
on different scales as provided by a filtering of the cluster density fluctuation 
field by a simple Gaussian smoothing filter. This function is shown in Fig. 
7 where we note that at the scale of the maximum of the power spectrum 
(a scale of ~ 100 h^ Mpc) the cluster density varies typically by about 
10%, while at the largest scales (400 /ij~ g Mpc) a fluctuation level of the 
order of 1% is indicated. This sets a high standard of requirements on the 
cluster survey, and while we are confident that the present sample contains 
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Figure 6: Power spectrum of the density distribution of the galaxy clusters 
in the REFLEX sample. For the analysis only 188 clusters in a box of 400 
Mpc side length was used (see Schuecker et al. in preparation). A value 
of Hq = 100 km s _1 Mpc -1 is assumed for the scaling. The solid line is a 
parametric fit to the data not considering a particular cosmological model. 

no systematic biases on the 10% level, further tests and simulations are 
needed to explore the reliability of the results on the largest scales. 

4 Conclusions 

It is obvious from the above arguments that studies of the large-scale struc- 
ture using clusters as probes require great care in the preparation of the 
sample of test objects. The first results of the REFLEX Survey presented 
here and some further tests that we have conducted already demonstrate 
the high quality of the present sample. This is the result of a very ho- 
mogeneous and highly controlled selection process used to construct the 
sample. The challenge of the selection work in this survey is to combine 
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Figure 7: Square root of the variance of the fluctuations in the cluster density 
distribution as a function of scale (obtained with Gaussian filtering) . A value 
of Hq = 100 km s" 1 Mpc -1 is assumed for the scaling. 

the complementary information from X-ray and optical wavelength in the 
most homogeneous way. Contrary to several earlier studies that we have 
conducted (Pierre et al. 1994, Romer et al. 1994, Ebeling et al. 1997, De 
Grandi 1996) in which the identification process was optimized to find as 
many clusters as possible by using all available sources of information, we 
have now achieved a highly complete cluster selection by just combining the 
ROSAT All-Sky Survey data and the COSMOS optical data base in a ho- 
mogeneous way, completely controlled by automated algorithms. Additional 
information is only used in the final identification but does not influence the 
selection. This is a very important achievement in this survey work. 

The data presented are still not fully complete in redshifts. But with 
data already obtained in January and September 1998 we can practically 
complete this data set (to 96%). An extended sample of REFLEX clusters 
down to a flux limit of 2 • 10" 12 erg s" 1 cm" 2 is already prepared and 
redshifts are available for more than 70% of the objects. This extended set 
of about 750 galaxy clusters will help very much to tighten the constraints 
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for the power spectrum and extend it to larger scales. It will further enable 
us to investigate the cluster correlation function - in particular the X-ray 
luminosity dependence of the clustering amplitude, which is an issue not yet 
resolved. Finally, a complementary ROSAT Survey cluster identification 
programme is being conducted in the Northern Sky in a collaboration of 
the Max-Planck Institut fur extraterrestrische Physik and J. Huchra, R. 
Giacconi, P. Rosati and B. McLean which will soon reach a similar depth 
and provide an all-sky view on the X-ray cluster distribution. 
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